Deciphering Inflation with Gravitational Waves: Cosmic Microwave Background 
Polarization vs. Direct Detection with Laser Interferometers 



o 
o 

cn 
> 

m 
O 

o 

:^ 

6 



X 



Tristan L. Smith^, Hiranya V. Peiris^*, Asantha Cooray'^ 
^California Institute of Technology, Mail Code 130-33, Pasadena, CA 91125 
^Kavli Institute for Cosmological Physics and Enrico Fermi Institute, University of Chicago, Chicago IL 60637 
Center for Cosmology, Department of Physics and Astronomy, University of California, Irvine, CA 92697 

(Dated: February 5, 2008) 

A detection of the primordial gravitational wave background is considered to be the "smoking-gun" 
evidence for inflation. While super-horizon waves are probed with cosmic microwave background 
(CMB) polarization, the relic background will be studied with laser interferometers. The long lever 
arm spanned by the two techniques improves constraints on the inflationary potential and validation 
of consistency relations expected under inflation. If gravitational waves with a tensor-to-scalar 
amplitude ratio greater than 0.01 are detected by the CMB, then a direct detection experiment with 
a sensitivity consistent with current concept studies should be pursued vigorously. If no primordial 
tensors are detected by the CMB, a direct detection experiment to understand the simplest form of 
inflation must have a sensitivity improved by two to three orders of magnitude over current plans. 



I. INTRODUCTION 

Recent high precision measurements of the cosmic mi- 
crowave background (CMB) anisotropy power spectrum 
01 have confirmed inflation j2j as the dominant paradigm 
to explain the origin of primordial fluctuations with a 
nearly scale-invariant spectrum. In addition to density 
perturbations, inflationary models predict a stochastic 
background of gravitational waves [j, 01 with the ampli- 
tude of the gravitational wave background given by the 
height of the inflaton potential when relevant modes exit 
the horizon. 

Inflationary gravitational waves with wavelengths com- 
parable to the horizon size are now being sought via 
ground-based CMB polarization experiments and 
eventually with a dedicated satellite generally named 
CMBPol 6]. In addition to the CMB effort, concept 
studies are now underway to investigate the possibility 
of directly detecting the relic background with a laser 
interferometer in space (Big Bang Observer 0; for the 
DECIGO proposal in Japan, see Ref. 0). The direct de- 
tection technique will be sensitive to modes with wave- 
lengths roughly an arm length of the interferometer. 
Based on the expected foreground confusion and tech- 
nological improvements, current concept studies aim for 
the frequency regime between 0.1 Hz to a few Hz. Since 
physical scales probed by the CMB and laser interferom- 
eters differ by ~17 orders of magnitude in wavelength, 
the large lever arm produced by combining the two tech- 
niques allows the inflaton potential to be pinned down 
better than any single method Hj, (lOj, (111, O • 

Here, we consider two separate studies on inflation us- 
ing the combined information from CMB and direct de- 
tection experiments. First, based on a Monte Carlo de- 
scription of the inflationary dynamics, we study the rel- 
ative abilities of the CMB and a direct detection method 
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to probe the inflaton potential in detail by making no as- 
sumptions on the power-law behavior Il0| or on a model 
shape for the potential 11] . We flrst consider inflationary 
models allowed by CMB data, making use of constraints 
on both the scalar and tensor spectrum, and address if 
laser interferometers can further improve the identifica- 
tion of potentials. For the CMB, we make use of the 
expected level of uncertainty with Planck 30|| , a possible 
detection with CMBPol, and the foreground-limit from 
CMBPol. For details on the potential detectability of the 
tensor amplitude with CMB polarization observations, 
we refer the reader to Ref. [2l|. For the direct detec- 
tion experiments, we make use of predictions related to 
BBO and DECIGO 0. Since our goal is to see how 
the simplest models of inflation can be constrained, our 
comparisons are for a general single-field potential in the 
slow-roll regime. 

In the second part of the paper, we drop all assump- 
tions related to single-field slow-roll inflation and study 
how well the CMB and the direct detection experiments 
can be combined to constrain the singic-fleld consistency 
relation between tensor spectral index and the ratio of 
tensor-to-scalar amplitudes for single-fleld inflationary 
models. Previous work on this possibility is found in 
Ref. [2^ where the test was limited to simply the in- 
formation on the tensor spectral index from CMB data 
alone. We show the extent to which a CMB-only analysis 
can be improved by adding direct-detection information 
as the latter allows a better determination of the tensor 
spectral index. 

This paper is organized as follows: In the next Section, 
we present the Monte Carlo modeling of the inflationary 
potential that satisfy both detections or limits at CMB 
and laser interferometer scales. In Section III, we discuss 
a fundamental test of single-field slow-roll inflation, in- 
volving the consistency relation between tensor-to-scalar 
ratio and the tensor spectral index. We conclude with a 
brief summary of our results in Section IV. 
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II. MONTE CARLO CALCULATIONAL 
METHOD 

To calculate observable spectra, we make use of a 
Monte Carlo technique to formulate the inflationary dy- 
namics through an infinite hierarchy of fiow equations 
involving the generalized "Hubble Slow Roll" (HSR) pa- 
rameters 

[illlllllll3- We link the Hubble parameter 
directly to the field (j) instead of time, H = H{(f)), un- 
der the assumption that is monotonic in time. The 
equation of motion for the background is 
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where the inflaton fleld evolves as = — mp[7J'((/))/47r. 
Here, an overdot corresponds to the time derivative and 
a prime denotes the derivative with respect to (j). The 
advantage of this formulation is that one can study the 
generic behavior of slow roll single fleld inflation without 
assuming a particular shape for the potential, except for 
the assumption of a single fleld. In terms of the HSR 
parameters ^Xh, the dynamics of inflation is described 
through: 
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Substituting equation lO in equation gives the infla- 
ton potential 
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The trajectories of the flow parameters are now gov- 
erned by a set of coupled first order differential equations. 
In practice, one has to truncate the infinite hierarchy at 
some finite order; in this paper we retain terms up to 10th 
order. Truncating the hierarchy of flow parameters at the 
term Xh means that ^^+^A^f = at all times as well. 
From Eq.El it also follows that = 

at all times. This simply describes a polynomial of order 
M + 1 in H{(l)) [13 with 
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H{cf) = Ho 
Further, from the deflnition of e(0). 



M+l' 



(5) 



when the coefficients Ai, with i > 1, are written in terms 
of the initial values of the HSR parameters as 
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where Ai = y/Aneo specifies the direction the field is 
rolling. These slow roll parameters require a prior as- 
sumption on the ranges of values taken. In the absence 
of any a priori theoretical knowledge, one can assume flat 
priors with the requirement that the potential satisfies 
the slow-roll condition; the latter is simply a statement 
about the smoothness of the potential. Since we are not 
attempting to make any statements about the measure 
of infiationary trajectories, but simply use the method 
as a Monte Carlo generator for potentials satisfying the 
slow roll conditions, such an assumption is justified. This 
work does not depend on the density of potentials as a 
function of tensor-to-scalar ratio, which is determined by 
the measure on the initial conditions. 

Without loss of generality, we can pick some fiducial 
physical scale that corresponds to 4>cmb^ which we take 
to be at kcMB = 0-002 Mpc~^. Then, with the above 
convention, (j) > 4>cmb corresponds to scales larger than 
^CMB (i-c. going further back in time), and (j) < 4>cmb 
corresponds to smaller scales. The physical wavenumber 
is associated with a value of <j) through 
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while the number of e-folds before the end of inflation, TV, 
comes from d4>/dN = mpi^/e/2y/Tr, with the convention 
that N increases as one goes further back in time. Here, 
we require that each potential generated by the Monte 
Carlo flows provide at least = 55 e-folds of inflation. 
For potentials where infiation ends through the break- 
down of slow-roll, the CMB observables are calculated 
at 55 e-folds before the end of infiation; for potentials 
corresponding to the hybrid case, the CMB observables 
are arbitrarily calculated at the 600th e-fold, assuming 
inflation ends at 655 e-folds through an orthogonal mech- 
anism. 

The standard observables are given in terms of the flow 
parameters to second order in the slow roll 0, |23| : 

Hs = 1 + 2t]~ Ae~2{l + Cy (9) 
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where ris, r, and rit are the tilt of the scalar spectrum, 
tensor-to-scalar ratio, and the tilt of the tensor spectrum, 
respectively. Additionally, we consider the running of the 
scalar tilt with 



dris 
din k 



dN 



[A + A{l + C) e] 



de_ 

dN 



(12) 



(3 - 5C) / dr] d€_\ (3 - C) 1 



dN 



'dN 



dNi ' 



3 




10" 
10" 

,E 10" 
E 

\ 10" 
> 10" 
10" 
10" 



-2.0 -1.5 



1.0 -0.5 




BBO sensitivitv floor 




- DECIGO sensitivity floor X lO^^^^^H 









-0.6 -0.4 -0.2 0.0 0.2 0.4 
* / mpi 



FIG. 1: The set of potentials satisfying 0.9 < < 1.0 from the Monte Carlo flow simulations: Left: constraints from Planck 
■with r = 0.02±0.01, dus/dlnk — O.OiO.Ol, BBO-standard, BBO-grand (factor 10 more sensitive than BBO-standard) Center. 
constraints from CMBPol "with optimistic foregrounds 21]: r = 0.001 ± 0.0003, dn^/dlnfc = 0.0 ± 0.005, BBO-standard, BBO- 
grand, and Right: the sensitivity limit due to foregrounds for CMBPol with r < 10^* and DECICO (factor 400 more sensitive 
than BBO-standard). Color coding: red and blue denote the CMB experiment without and with the dUs/dlnk constraint, 
respectively. Creen: BBO-standard. Black: BBO-grand {Left, Center) and DECICO (Right). The direct detection constraints 
are applied to the tensor amplitude and tilt at 1 Hz, following the procedure described in the text. If a particular color does not 
appear in a plot, it has been overwritten by the next tightest constraint, and the latter is therefore not helpful in constraining 
the potential. The meaning of the color coding is further clarified in the text. Here (p — corresponds to CMB scales while 
curves end at 4> < corresponding to a frequency of 1 Hz probed by direct detection methods. 



where de/dN = 2e(ry-e), dr]/dN = -€!]+£_, and d£_/dN = 
^(?7 — 2e)-f ^\h- Here and above, C = 4(ln 2 + 7) — 5 and 
7 — 0.5772 is the Euler-Mascheroni Constant, 77 = ^A^f, 
and ^ = ^Aff. 

For direct detection experiments, we take a fiducial fre- 
quency of / = 1 Hz with fcdir = 6.47 x IQ^^ Mpc"! for 
observations. With ln(fcdir/fccMB) = 40.3 the large lever 
arm is expected to improve constraints on the inflation- 
ary model 1^ . The direct detection observables are cal- 
culated from the potentials by finding the (j) correspond- 
ing to ln(A:dir/fccMB) = 40.3. In these experiments, at 1 
Hz, the signal-to-noise ratio for a detection of the grav- 
itational wave background is SNR = X [VIq^ / ^ 
where X is '--^ 0.25, 2.5, and 100 for concept study de- 
signs involving a standard BBO, an optimistic version of 
BBO, and DECIGO. The last two possibilities improve 
sensitivities through multiple detector correlations. The 
uncertainty of the tensor spectral index at 1 Hz is taken 
to be ariT ~ 6/SNR [3. 

Fig- n shows a set of potentials from the Monte Carlo 
flow simulations that satisfy levels of tensor and scalar 
modes reachable by Planck, a CMBPol-style polarization 
satellite designed to probe primordial gravity waves [IJ , 
and the foreground limit of this satellite. We show all 
potentials with 0.9 < < 1.0, though in practice the 
CMB measurements will also yield strong constraints on 
this parameter. In Fig. 1 we also show the improvement 
in constraints on the potential using information from di- 
rect detection experiments with uncertainties in the ten- 
sor amplitude and tilt calculated following results of the 
analysis in Ref. [T^ . 

In Fig. 1, the potentials generated by the Monte-Carlo 
process are color-coded in the following way. First, all po- 
tentials with 0.9 < ris < 1.0 which are not ruled out by 
applying the constraints on r from the CMB experiment 
are shown in red. Then the potentials which are not ruled 



out by further applying the dn^/d In k constraint from the 
CMB experiment are over-plotted in blue. Next, the po- 
tentials which are not eliminated by adding the constraint 
from BBO-standard on the tensor amplitude and tilt to 
the CMB experiment are over-plotted in green. Finally, 
the potentials which are not eliminated by adding the 
constraints from either BBO-grand (left, centre) or DE- 
CIGO (right) to the CMB experiment are over-plotted in 
black. At each stage, the elimination of a larger fraction 
of potentials indicates the usefulness of the extra con- 
straints. If a particular color does not appear in a given 
panel, it means that all potentials have been overwritten 
by the next tightest constraint, and therefore the next 
tightest constraint does not aid in constraining the po- 
tential. Therefore, while the combination of Planck and 
BBO-standard, and especially BBO-grand, leads to an 
improvement in constraints on the potential, these di- 
rect detection experiments do not compete at all with 
CMBPol either in the case of a detection at r = 0.001 or 
the limit of r < 10~^ at CMB scales. A more sensitive 
experiment such as DECIGO probing Hqt^H^ > 10~^° 
competes well with CMBPol if r < 0.01 and is both de- 
sirable and useful to understand inflation. 



III. A FUNDAMENTAL TEST OF SLOW-ROLL 
INFLATION 

To further understand the usefulness of a direct detec- 
tion experiment, we also study how well the inflationary 
single- field consistency condition can be tested |2a|. as 
any departure can capture important physics |27l| . We 
define TZ = —r/8nt so that TZ = I corresponds to the 
consistency relation. The uncertainty in determining TZ 
(given 7?, = 1 as the fiducial value) can be written in 
terms of uncertainty in measurements of r and nt assum- 
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FIG. 2: Left: The mapping between tensor-to-scalar ratio r at CMB scales and Qcwh^ at 1 Hz for a direct detection. The 
gray shaded region shows the uncertainty implied with = 0.95 ±0.1 [2^ by keeping terms up to second order (running of the 
tensor spectral index) in the slow-roll power-law expansion. We include the line 'Unphysical for inflation' to indicate the region 
above which rit > 0. The three horizontal curves are the 2<j detection amplitudes for BBO-standard, BBO-grand, and DECIGO 
in solid, dashed, and dotted lines respectively. In these three experiments, at 1 Hz, the signal-to-noise ratio for detecting the 
gravitational wave background is SNR — X(riGw/10~^*), with values for X shown in the panel. Right: The la uncertainty 
in the single-field consistency relation TZ = —r/8nt. The thin lines follow laser interferometers in the right panel, showing the 
error expected by combining direct detection measurements of rit with CMBPol measurements of the tensor-to-scalar ratio, r. 
The two thick lines indicate errors in TZ when rit is determined from CMB alone. The thick solid curve corresponds to the 
expected accuracy of ESA's Planck satellite, the thick dashed curve corresponds to CMBPol. The shaded region indicates a 
50% determination of the consistency relation {TZ — 1.0±0.5). For direct-detection observations, sensitivity is degraded at large 
r because of an increase in the uncertainty of the importance of the running of the tensor spectral index from direct detection 
scales to CMB scales; at small r, the accuracy with which TZ can be determined is dominated by the error in measuring rit with 
CMB '^^l and laser interferometers where the latter is cr„j oc 1/r. 



ing these are uncorrelated, 
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A 10% determination of TZ with the expected value 
of unity requires measuring rit with an uncertainty of 
cr„, ~ 0.0125r. Using tables of Ref. 21], CMBPol with 
a 10% foreground contamination and r = 0.01 gives 
TZ = 1.0 ± 80. The obstacle is the inabihty of CMB 
polarization observations to measure rit precisely, since 
polarization anisotropics probe a limited range in the un- 
derlying tensor spectrum modes, and the range probed 
is also contaminated by cosmic shear p^ . 

When combining with a direct detection experiment, 
however, the situation improves significantly. We calcu- 
late the la error expected in TZ at CMB scales based on a 
determination of nt from direct detection scales and a de- 
termination of r by CMBPol. In Ref. [l3| it was shown 
that a laser interferometer (whose sensitivity peaks at 
~ 1 Hz) would be able to determine the spectral tilt of a 
gravitational wave background, with an amplitude r^cWj 
to an accuracy of 
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where the various values of X corresponding BBO- 
standard, BBO-grand and DECIGO are shown in Fig. 2 



and h is the Hubble parameter today in units of 
100 kms-iMpc"^ 

While we take our fiducial model to be single-field 
slow -roll inflation (i.e. TZ = 1.0), we make no assump- 
tions about single-fleld slow-roll inflation when we re- 
late f^GW^^ to the tensor-to-scalar ratio. First we relate 
flowh^ to the primordial spectrum as in Ref. jlli |. 



(15) 



where Pt{k) is the primordial power spectrum of infla- 
tionary gravitational waves, and Aqw — 2.74 x 10~^. 
The factor, Aq-^/v, takes into account how the gravita- 
tional waves have evolved after re-entering the horizon. 
The primordial inflationary gravitational wave power 
spectrum can be approximated as a power law with a 
running spectral index, 



Pt{k) « Pt(fco) ( ^ 



^ \ nt + ^Qt ln(fc/fco) 



(16) 



Using the expressions in Eqs. H9I11|I to first order in the 
slow roll parameters along with the expression for the 
running of the tensor spectral index, at, 
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we are able to express rit and at in terms of r and rig 
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where Ps{ko) is the amplitude of the scalar perturbations 
at some pivot wavenumber ko and all of the spectral 
quantities are also measured at this wavenumber. For 
this analysis we take ru — 0.95 and Ps(^o) 2.21 x 10^^ 
at fc = 0.002 Mpc"^ (if. 

Since we must connect a measurement of nt at BBO 
scales to the rit at CMB scales in order to determine 
TZ, we must posit some scale-dependent relation between 
these two measurements. We emphasize that we can- 
not use the flow equation approach of the previous sec- 
tion here; those equations assume single-field inflation 
and therefore implicitly embody the consistency relation, 
whereas here we are attempting to test the single-field 
assumption. Instead, we assume that the gravitational 
wave spectrum is close to scale invariant with an un- 
known, but higher order, running so that nf^'^ « np^^. 
Since the running of nt is unlikely to be determined 
by either the CMB or direct detection we include an 
additional uncertainty due to the unknown running, ^ 
(1/2) ln(/cdir/fccMB)"?, as in Eq. ((T7I) . This leads to a 
decrease in the sensitivity of direct detection measure- 
ments of TZ at 'large' values of r. In particular, the error 
on our determination of nt is given by the expression, 
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Our results are summarized in Fig. 2. The left panel 
shows the mapping between tensor-to-scalar ratio r at 
CMB scales and JIgw^^ at 1 Hz for a direct detection. 
After this paper appeared in preprint, another study p9| 
was submitted, where a criticism was made of a region on 
this panel, describing it as unphysical. In fact this region 
is only unphysical if inflation is being assumed as the gen- 
erating mechanism for tensor modes, and one cannot do 
that if one is attempting, as we are, to test that assump- 
tion in the first place. Therefore we continue to show 
this region, while indicating that it is unphysical under 
the assumption of inflation. As shown in the right panel, 
BBO improves relative to CMBPol alone by a factor of ^ 
a few in the uncertainty of 7?. if r ~ 0.1. This is unlikely 
to be useful given that current observations already limit 
r to be below 0.3. A determination of 7?. as 1.0 ± 0.5 is 
achievable when 10"^ > r > 10"^ with DECIGO, while 
if r ~ 10~^, it is unlikely that even DECIGO would pro- 
vide a determination of 7?, to a reasonable accuracy. In 
general, either version of BBO is unlikely to be competi- 
tive with CMBPol, and the sensitivity level of DECIGO 
must be considered as the experimental target goal to 
pursue a direct detection at ^ 1 Hz. Further limits on 



r from the CMB will tighten this conclusion and could 
only lead to a further increase in required sensitivity un- 
less a detection is made with the CMB at r > 0.01. Our 
conclusions are independent of the choice of a flducial 
frequency between 0.1 Hz and few Hz, but could be sub- 
jected to the highly uncertain impact of foregrounds at 
direct detection frequencies Additional physics be- 
tween the CMB and 1 Hz scales only strengthen our 
conclusions on the required sensitivity for a laser interfer- 
ometer as these exotic models generally lower the tensor 
amplitude further. 



IV. SUMMARY 

To summarize, we have considered the relative 
strengths of CMB polarization observations and direct 
detection laser interferometers in constraining the infla- 
ton potential. For single-fleld slow-roll inflation models, 
without relying on any particular shape for the poten- 
tial, we find that direct detection experiments with sen- 
sitivities around BBO can improve constraints on infia- 
tionary models relative to Planck. However, when com- 
bined with CMBPol, these direct detection sensitivities 
are unlikely to be competitive. While we have not con- 
sidered exotic models, and a case can certainly be made 
for a low-sensitivity direct detection experiment based on 
non-standard descriptions for inflation including models 
of bubble nucleation and pre-Big Bang descriptions 
pif. it is also important to understand first how these 
experiments test the simplest forms of inflation. 

In this context, we also discuss a determination of the 
single-field slow-roll consistency relation, which is a way 
to establish an underlying model within the inflation- 
ary paradigm to probe physics at the earliest times of 
the Universe. In general, we find an experiment like DE- 
CIGO, with sensitivity level of ^cwh'^ > 10^^" to be the 
preferred option; however, it is unlikely that the consis- 
tency relation will be determined to the accuracy needed 
to see loop corrections 27], unless the tensor-to-scalar 
ratio is greater than 0.05 and loop corrections are at the 
level of 10% or more. 
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